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ABSTRACT 

The problem of transient heat transfer through a monatomic collisionless 
ideal gas enclosed between parallel walls with wall accommodation coefficients 
of unity is analyzed for the case of a step change in one wall temperature. 

The results indicate that steady state heat transfer condition is ap- 
proached in an oscillatory manner requiring about 20 to 30 crossings of the 
atoms across the channel at the mean thermal velocity based on the initial 
wall temperature, C wi = (2R g T^) 1 / 2 , to achieve steady state. 

The solution is obtained in terms of an arbitrary quantity Q transported 
across a unit elemental area some distance X from the channel wall. Results 
for other transported quantities than energy i.e. density, temperature, can 
be obtained from the present formulation. The results can be readily extended 
to an arbitrary change in wall temperature as a function of time. 

INTRODUCTION 

The interest in rarefied gases has been increasing largely due to the 
possibility of flight at high altitudes and high speeds. This is because 
the mean free paths for molecular collisions increases with increasing altitudes. 
For instance (ref. l) at an altitude of 62.5 miles the mean free path 
is 0.1 ft. The mean free path increases to 100 feet at an altitude of 
120 miles. Thus, convective heat transfer in a earth satelite open to the 
atmosphere would be through a rarefied gas. 
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There has been very little consideration given to transient heat trans- 
fer between surfaces enclosing a collisionless gas, i.e. where intermole cular 
collision of the molecules are sufficiently rare so that they can be neglected. 
The problem of transient heat transfer from a single plate to a collisionless 
gas has previously been treated by Yang and Lees (ref. 2) for a step increase 
in wall temperature. They also included the effects due to a step increase 
in wall velocity in a direction in the wall plane. 

The present analysis considers the effect of a step change in one wall 
temperature on heat transfer when the collisionless gas is enclosed between 
two walls. In this case the Yang and Lees solution would apply for early 
times before the effect due to the second wall would be felt. In the analysis 
we will determine how long and in what manner the final steady state condition 
is achieved. If the time for the transient to be completed is a significant 
part of the time used in a heat transfer situation then this transient effect 
would have to be included in the analysis. 

The present analysis can be readily extended to the case of an arbitrary 
change in wall temperature with time. In this case if the rate of change in 
the wall boundary temperature is small compared to the time needed for steady 
state conditions to be attained then the transient effects also would have 
to be accounted for in calculating the heat transfer. This assumes that the 
convective heat transfer is significant compared to the radiant heat trans- 
fer. The radiant heat transfer transient for this same problem takes place 
almost instantaneously so that the transient can be neglected for that case. 

The present analysis also serves as a limiting solution to the more 
general and more difficult problem of transient heat transfer in the case 
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of a rarefied gas where the intermolecular collisions are sufficiently 
numerous so that this effect must he included in the analysis. 

The present results show the time needed for steady state conditions 
to he achieved and the behavior of the heat flux as it changes in time from 
its initial value to its final steady state value. These results are shown 
for several different values of the ratio of final to initial wall temperatures. 
Results for hoth step increases and step decreases in wall temperature are 
shown. It was assumed that the wall accommodation coefficient was unity and 
the gas was ideal and monatomic. 

The solution is obtained in general form so that various other quantities 
in the gas such as temperature and density could he obtained as was done in 
Ref. 2. These results were not carried out because it was felt that they 
would not be significant in a collisionless gas. 

NOMENCLATURE 

C average thermal velocity (2RgT)' L / 2 

D distance between walls 

M mass of a atom 

mass flux leaving the wall, (pV 2 ) w = (p^^^it 1 / 2 ) 

Q molecular quantity 

q total energy flux 

R 2 ratio of wall temperatures, (T /T .) 

c 7 W,0' W,l' 

Rg gas constant 

r see fig. 1 

S see fig. 1 

T temperature 

u inverse dimensionless velocity, l/v 
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V molecular velocity 

v dimensionless velocity, V/C^ 

X coordinate normal to wall 

6 dimensionless quantity, rC^/X 

9* dimensionless quantity, tC^/(D - X) 

| ratio of densities, p/p± 

p gas density 

t time 

f dimensionless time, tC^/D 

() refers to averaged quantity 

Subscripts: 

+,- to or from wall wo 

i initial conditions before transient occurs 

s steady state 

wo,wl refers to wall wo or wall wl 

Superscript: 

refers to quantity leaving wall at time t' 

ANALYSIS 

The problem treated here is the heat transfer through a collisionless 
gas enclosed between two infinite parallel plates when the temperature of one 
plate is changed in a stepwise manner. Of prime interest is the net heat 
transfer between the gas and walls as a function of time from the initiation 
of the transient until steady state heat transfer conditions are achieved. 

The problem of transient heat transfer from a single plate to a col- 
lisionless gas has previously been treated by Yang and Lees (ref. 2). They 
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also included the effects due to a step increase in wall velocity in a 
direction in its own plane. The results of reference 2 for one wall agree 
with the present result for the wall that has its temperature suddenly 
changed at early time before the effects due to the opposite wall are felt. 

The analytical model is shown in figure 1. Two plates enclose a col- 
lisionless monatomic ideal gas, i.e. , the mean free path to collision between 
molecules in the gas is large compared to the distance across the channel. 

The plates are assumed to have an accommodation coefficient of 1, that is 
all the incident molecules come to equilibrium with the wall during collision. 
Initially the entire system is at temperature T^. At time t = 0 there 
is a step change in the wall wo temperature from T^ to T wo . The upper 
wall temperature remains fixed at T wl = T wi' 

The analysis of the present problem proceeds as follows: The amount of 
some molecular quantity Q(V) transported across an elemental area at X at 
time t that originated at the lower surface elemental area dA wo at time 
t' (fig. l) can be derived from reference 3 equation (Al) to be 



where the mass flux leaving the wall, m^ 0 is given by ( Pw 0 C W0 ) . 

Taking a ring element dA' = 2itr dr dr with S 2 = X 2 + r 2 equation (l) 


can be rewritten as 
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The molecules that reach the elemental area dA at X will he of two kinds. 
One group will have left the wall before the change in the wall temperature. 

The second group reaching the elemental area will have left after the transient 
is initiated. Considering the first group, part of these molecules that 
reach dA at X at time t that have velocities less than X/t for all 
values of path length S will have left the surface before the change in 
wall temperature . This term is given by part I in equation (3) where 6 = 
ty/X. For the molecules with velocities greater than X/t, only those mole- 
cules with path length S greater than xV will have left the surface before 
the change in wall temperature. This term is given by part II in equation (3). 

Now consider the second group of molecules arriving at dA at time x. 
This group has left the wall after the change in wall temperature and will 
consist of molecules with velocities greater than X/x and path length S 
less than xV. This term is given by part III in equation (3). 



Ill 
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Transforming variables, S 


P ) /u, v = l/u using the Jacobean gives 


dS dv 


1 

u 

0 


f _ 



df du 


dp du 


u" 


( 4 ) 


Combining equations (3) and (4) we obtain 



^ 0 Q(exp - (Ru)" 2 ) 
(Ru ( 6 - e*)) 3 


dO 1 du 


By changing the order of integration this can be written as 


(5) 


2^/2 

Pi C i 


(p(V 2 Q))+ = 4 




Q(exp - u~ 2 )du d9' 
u 3 (e - e') 5 


+ 4 



i; 0 ^ exp " ( Ru )' 2 ) 
rV (6 - e -) 3 


du d0* 


( 6 ) 


Similarly the net Q transported from the opposite wall wl at X of D 
across the differential area 4A at X is given by 
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2* 


1/2 


p'=o n u=e*-e*' 


--- c ~ ( p ( v ? Q >) = 4 
Pi c i 2 “ 


+ 4 



0*' =-co u=0 

e*'=e*r'u=e*-e*> tt A/ „-2, 


Q(exp - u~ 2 )du d0*' 
( 0* - 0* ' ) 3 u 3 


£^-|Q(exp - u )du d0*' 
(0* _ 0*') 3 u 3 


(7) 


F 0*'=O kJ u=0 

Then the net quantity Q transported across the elemental area at X is 
given by 


Pi c i 

( p(VgQ) ) =7^2 


p(V 2 Q) + + p(V 2 Q). 


( 8 ) 


Calculation of mass flux across the channel . - Assuming Q = 1 we can 
obtain the mass flux transported across an elemental area from equations (6) 
to (8). Carrying out part of the integration in equations (6) and (7) gives 


< p < V > + - 1 ^ 


(p(v 2»- ■ 1 - *‘ 1/e * 2 + 2 



0'=0 e « 

5xo 


WO 



i 1 

exp - - 

_R( 0 - 0' ) 

0 



0'=O 


0.-0. ^ 


r( 0 - 0’) 3 


d0' 


^ xp (0. - 0. 


■) 2 J 


(0* - 0*')' 


(9a) 

d0*' 

(9b) 


'0*'=O 

Since at a wall the incident mass flux is equal to the reflected mass flux 

(p<V 2 » +u = (p(V 2 ».„ (10) 

Then after the wall transient has occurred, t > 0, at walls wO and wl we 
have 


( P< v 2» t ,vi < e < Vh ^ 

(c(V )wi ° (p(VVi = 6¥l 


(11a) 
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(p(v 2 »- (p<v 2 ))+ ,WO 

( p(^ 2 ) )vi ( p(Vg) ^ wo 


(lib) 


Combining equations (9a) with (lla) and (9b) with (lib) we obtain a pair of 
simultaneous integral equations 


6vl = 1 " 


.f-2 


W* * 1 ' 



exp |j^R(^ - ^ , )}"^jdt' (12a) 


)T(t - f)"^] 


exp -(t - P)~ d dp (12b) 


These are solved by an iterative numerical integration procedure on a high 
speed electronic computer and the results are shown in figure 2. Also shown 
in figure 2 are the limiting results for very large times, P*x>. For this 
case, equations (12a) and (12b) can be shown to give 


^vl,P>co = 1 + R ( 15a ) 

2R 

R two,P°° = 1 + R ( 13b ) 

This limiting result for large time agrees with the steady state collisionless 
solution given in reference 4. For very short times, P+0 equations (12a) 
and (12b) become 


^Wl,P+0 : 1 

R ^vo,P+0 = 


(14a) 

(14b) 


This limiting result for small times agree with the results of reference 2. 


Energy flux from the channel wall to the gas . - For a monatomic gas the 
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energy transported across the channel is due to translational motion only 
(ref. 2). The total energy being transferred across the channel is therefore 
given by 


* ■ •(••©') 

Then in equations (6) to (8) letting Q = V 2 /2 = cS/2u 2 we obtain 


2jc 

P 


. U . exp I. iV ^ ^ 

± c[ \ e* 2 ) \ e 2 / 2 e* 2 2e 2 



0* —Q 


R| wo expf- 


0'=O 


(9 

9*'=©* Iwi exp 


, ^( e 2 e> ) fi +[ r(s - e'Q^de* 


( 0 * 


r 0*»=O 


(e* _ e*>)3 


fear, 


-§ 


(15) 


+ (0* _ 0*')- 2 d0* (16) 


To find the total heat flux leaving the walls we write 0 and 0* in 
terms of X and f as given in the nomenclature and take the limits for 
X = 0 and D. This gives at wall wo the result 


,S ^j o ' r3s - w ' 1 + 2 ( L + i?) 




wfr 4 ( * ; !2 '3 If . 


/p=0 


(f - f ’ ) 


m3 


(17a) 


Similarly the result at wall wl is given by 
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l - ! wl (t) 



Ri wo (t'){l + [R(t - t'j]" 2 ) 


Afr»=0 


(i]r - ilr * ) 1 


exp 


- (r(# - t’)) 


-2 


dp (17b) 


These results are shown in figures 3(a) and (b). 

The limiting solution for large times, poo is the fully developed 
collisionless result and agrees with the value given in reference 6. 


^/2 C 


Pi 




= R(R - 1) 


J t=°°,X=0,X=D 


(18) 


For short times, pO we obtain the limiting results for early time for 
wall wo which agree with reference 2 



po,x=o 



(19) 


while at wall wl for early times the heat transfer from the wall to the gas 
is zero 


1/2 

ix ' 


Pi 


C 3 

°i 


= 0 


'PQ,X =: D 

These results are shown in figure 3. 


RESULTS AND CONCLUSIONS 


( 20 ) 


The results for the mass flux leaving the walls are shown in figure 2. 
Initially the mass flux leaving walls wo and wl remain constant after a 
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step increase in wall wo temperature. Then the mass flux leaving the wall 
wl first starts to change towards the steady state result. Somewhat later 
the wall wo also begins to change towards the steady state result. The 
results approach the final solution in an oscillatory manner. 

The mass fluxes leaving the walls do not change for a short time after 
the step change in wall temperature because although the wall wo temper- 
ature changes at t = 0 the incident mass flux on that wall remains constant 
and so the mass flux leaving the wall remains constant from continuity. 

However in the case of the step rise in wall temperature the molecules leaving 
the heated wall wo will travel across to wall wl Taster because of their 
higher thermal velocity causing the incident mass flux on wall wl to 
increase as shown in figure 2. This is because not only are the faster 
moving molecules from the hot wall wo now arriving but the slower moving 
molecules that have left wall wo earlier also continue to arrive at wall 
wl. Then since the mass flux incident and therefore also leaving the unheated 
wall will now be greater, the mass flux incident on the heated wall will 
become greater and the mass flux leaving the heated wall will now start to 
rise towards the steady state result. Similar arguments can be used in 
explaining how the cooled wall temperature case approach the equilbrium mass 
flux. 

The results for the heat flux from the wall to the gas is shown in 

figure 3(a) for the case where the wall wo temperature is raised. It can 

be seen that immediately after the step change in wall temperature there 
is a heat transfer from the heated surface wo to the gas. The other surface 

wl with its temperature unchanged does not absorb any heat until the molecules 
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that have left the heated surface reach surface wl. Both walls approach the 
final steady state value asymptotically and are within 1 percent of the value 
hy about dimensionless time f ~ 10. 

The results for a wall temperature decrease is shown in figure 3(b). 

Heat is absorbed into the walls at the cooled surface wo immediately. The 
other surface wl with its temperature unchanged does not transmit any net 
heat to the gas until the molecules leaving the cooled wall reach it. Both 
walls then approach the steady state result in a oscillatory manner taking 
somewhat longer to reach the fully developed result than in the heating case 
i.e. for T wo /T w ^ = l/l6 1 percent of fully developed heat transfer is 
reached in dimensionless time f ~ 30. 

The present formulation of the problem allows the calculations of various 
properties in the body of the gas such as density and temperature by choosing 
appropriate properties of Q in equations (6) to (8). This was done in 
reference 2 for the heat transfer from a single plate. However it was felt 
for the collisionless gas these properties would not be useful. 

To illustrate the length of time needed for the step transient to reach 

final equilibrium we can run out a sample calculation. Taking time to reach 

1 percent of steady state in seconds, temperature in degrees rankine and 

distance in feet, from the definition of f we can write time t as 

s s 


Id Id 

S S 



Taking the molecular weight, m.w. , of air equal to 29, the distance between the 
plates D as 2 feet and the initial wall temperature T • as 100° R then 



14 


if the dimensionless time # g needed to obtain fully developed conditions 
is set equal to 30 we find from (eq. 21) a time t ~ 0.1 seconds as about the 
time needed to achieve thermal equilibrium. This implies that for heat trans- 
fer changes on the order ®f seconds the transient collisionless heat transfer 
behavior can be significant. 

The times for the steady state heat transfer condition to be achieved 
could be significantly different if the mean free path to collision in the 
gas was smaller so that intermole cular collisions would have to be accounted 
for and also if the accommodation coefficient of the surface was not equal 1. 

The present results can be extended to arbitrary variations of wall tem- 
perature as a function of time by considering the dimensionless temperature ratio 
2 

R as a function of time in the analysis. 
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Dimensionless wall heat flux, q w /q w 



Dimensionless time, ^ = xCj/D 


(a) Ratios 2 and 16. 

Figure 3. - Dimensionless wall heat flux for different wall temperature ratios. 
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(b) Ratios 1/2 and 1/16. 


Figure 3. - Concluded. Dimensionless wall heat flux for different wall temperature ratios. 


